Metagenomic Analysis of Viromes of Aedes Mosquitoes across India.

Metagenomic analysis of Aedes aegypti and Ae. albopictus mosquitoes from diverse geographical regions of India revealed the presence of several insect viruses of human interest. Most abundant reads found in Ae. aegypti mosquitoes were of Phasi Charoen-like virus (PCLV), Choristoneura fumiferana granulovirus (CfGV), Cell fusing agent virus (CFAV), and Wenzhou sobemo-like virus 4 (WSLV4), whereas WSLV4 and CfGV constituted the highest percentage of reads in Ae. albopictus viromes. Other reads that were of low percentage included Hubei mosquito virus 2 (HMV2), Porcine astrovirus 4 (PAstV4), and Wild Boar astrovirus (WBAstV). PCLV and CFAV, which were found to be abundant in Ae. aegypti viromes were absent in Ae. albopictus viromes. Among the viromes analyzed, Ae. aegypti sampled from Pune showed the highest percentage (79.82%) of viral reads, while Ae. aegypti mosquitoes sampled from Dibrugarh showed the lowest percentage (3.47%). Shamonda orthobunyavirus (SHAV), African swine fever virus (ASFV), Aroa virus (AROAV), and Ilheus virus (ILHV), having the potential to infect vertebrates, including humans, were also detected in both mosquito species, albeit with low read numbers. Reads of gemykibivirus, avian retrovirus, bacteriophages, herpesviruses, and viruses infecting protozoans, algae, etc., were also detected in the mosquitoes. A high percentage of reads in the Ae. albopictus mosquito samples belonged to unclassified viruses and warrant further investigation. The data generated in the present work may not only lead to studies to explain the influence of these viruses on the replication and transmission of viruses of clinical importance but also to find applications as biocontrol agents against pathogenic viruses.

Circulation of a single hepatitis A virus genotype IIIA with two distinct clusters in different states of India.

With the changing hepatitis A epidemiology in India, focal viral outbreaks are being reported from different parts of the country. This study presents Hepatitis A Virus (HAV) strain characterization (period 2009-2020) from 18 states of India. For that, blood and stool samples (n â€‹= â€‹280) were screened for HAV RNA and sequences for 5'non-coding and VP3 regions were generated from positive samples (n â€‹= â€‹68). Presence of a single IIIA genotype in all samples indicated IIIA being the only HAV genotype currently circulating in India. Interestingly, it was evident that these strains form two distinct groups suggesting independent evolution of these two clusters.

Abundance of Phasi-Charoen-like virus in Aedes aegypti mosquito populations in different states of India.

Mosquitoes are known to harbor a large number of insect specific viruses (ISV) in addition to viruses of public health importance. These ISVs are highly species specific and are non-pathogenic to humans or domestic animals. However, there is a potential threat of these ISVs evolving into human pathogens by genome alterations. Some ISVs are known to modulate replication of pathogenic viruses by altering the susceptibility of vector mosquitoes to pathogenic viruses, thereby either inhibiting or enhancing transmission of the latter. In the present study, we report predominance of Phasi Charoen-like virus (PCLV, Family: Phenuviridae) contributing to >60% of the total reads in Aedes aegypti mosquitoes collected from Pune district of Maharashtra state using next generation sequencing based metagenomic analysis of viromes. Similar results were also obtained with mosquitoes from Assam, Tamil Nadu and Karnataka states of India. Comparison of Pune mosquito sequences with PCLV Rio (Brazil) isolate showed 98.90%, 99.027% and 98.88% homologies in the S, M and L segments respectively indicating less genetic heterogeneity of PCLV. The study also demonstrated occurrence of transovarial transmission as seen by detection of PCLV in eggs, larvae, pupae and male mosquitoes. Ae. aegypti mosquitoes collected from Pune also showed a large number of reads for viruses belonging to Baculoviridae, Rhabdoviridae, Genomoviridae and Bunyaviridae families. The role of PCLV in the replication of dengue and chikungunya virus is yet not clear. It warrants further studies to know the significance of PCLV and other ISVs on the replication and transmission of Ae. aegypti borne pathogenic viruses, especially in the absence of prophylactics or therapeutics.

Phylogenetic classification of the whole-genome sequences of SARS-CoV-2 from India & evolutionary trends.

BACKGROUND & OBJECTIVES: Several phylogenetic classification systems have been devised to trace the viral lineages of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, inconsistency in the nomenclature limits uniformity in its epidemiological understanding. This study provides an integration of existing classifications and describes evolutionary trends of the SARS-CoV-2 strains circulating in India. METHODS: The whole genomes of 330 SARS-CoV-2 samples were sequenced using next-generation sequencing (NGS). Phylogenetic and sequence analysis of a total of 3014 Indian SARS-CoV-2 sequences from 20 different States/Union Territories (January to September 2020) from the Global Initiative on Sharing All Influenza Data (GISAID) database was performed to observe the clustering of Nextstrain and Phylogenetic Assignment of Named Global Outbreak LINeages (Pangolin) lineages with the GISAID clades. The identification of mutational sites under selection pressure was performed using Mixed Effects Model of Evolution and Single-Likelihood Ancestor Counting methods available in the Datamonkey server. RESULTS: Temporal data of the Indian SARS-CoV-2 genomes revealed that except for Uttarakhand, West Bengal and Haryana that showed the circulation of GISAID clade O even after July 2020, the rest of the States showed a complete switch to GR/GH clades. Pangolin lineages B.1.1.8 and B.1.113 identified within GR and GH clades, respectively, were noted to be indigenous evolutions. Sites identified to be under positive selection pressure within these clades were found to occur majorly in the non-structural proteins coded by ORF1a and ORF1b. INTERPRETATION & CONCLUSIONS: This study interpreted the geographical and temporal dominance of SARS-CoV-2 strains in India over a period of nine months based on the GISAID classification. An integration of the GISAID, Nextstrain and Pangolin classifications is also provided. The emergence of new lineages B.1.1.8 and B.1.113 was indicative of host-specific evolution of the SARS-CoV-2 strains in India. The hotspot mutations such as those driven by positive selection need to be further characterized.

The expression patterns of immune response genes in the Peripheral Blood Mononuclear cells of pregnant women presenting with subclinical or clinical HEV infection are different and trimester-dependent: A whole transcriptome analysis.

Hepatitis E is an enteric disease highly prevalent in the developing countries. The basis for high mortality among pregnant hepatitis E patients remains unclear. Importantly, a large proportion of infected pregnant women present with subclinical infection as well. In order to understand the possible mechanisms influencing clinical presentation of hepatitis E in pregnant women, we explored a system biology approach. For this, PBMCs from various categories were subjected to RNAseq analysis. These included non-pregnant (NPR, acute and convalescent phases) and pregnant (PR, 2nd and 3rd trimesters, acute phase and subclinical HEV infections) patients and corresponding healthy controls. The current study deals with immune response genes. In contrast to exclusive up-regulation of nonspecific, early immune response transcripts in the NPR patients, the PR patients exhibited broader and heightened expression of genes associated with innate as well as adaptive T and B cell responses. The study identified for the first time (1) inverse relationship of immunoglobulin (Ig) genes overexpression and (2) association of differential expression of S100 series genes with disease presentation. The data suggests possible involvement of TLR4 and NOD1 in pregnant patients and alpha defensins in all patient categories suggesting a role in protection. Induction of IFNγ gene was not detected during the acute phase irrespective of pregnancy. Association of response to vitamin D, transcripts related to NK/NKT and regulatory T cells during subclinical infection are noteworthy. The data obtained here could be correlated with several studies reported earlier in hepatitis E patients suggesting utility of PBMCs as an alternate specimen. The extensive, informative data provided here for the first time should form basis for future studies that will help in understanding pathogenesis of fulminant hepatitis E.

Detection of human parvovirus 4 DNA in the patients with acute encephalitis syndrome during seasonal outbreaks of the disease in Gorakhpur, India.

Seasonal outbreaks of acute encephalitis syndrome (AES) at Gorakhpur, India have been recognized since 2006. So far, the causative agent has not been identified. Use of next generation sequencing identified human parvovirus 4 (HPARV4) sequences in a CSF/plasma pool. These sequences showed highest identity with sequences earlier identified in similar patients from south India. Real-time PCR detected HPARV4 DNA in 20/78 (25.6%) CSF and 6/31 (19.3%) plasma of AES patients. Phylogenetic analysis classified three almost complete genomes and 24 partial NS1 sequences as genotype 2A. The observed association of HPARV4 with AES needs further evaluation. ELISAs for the detection of IgM and IgG antibodies against scrub typhus (Orientia tsutsugamushi, OT) showed ∼70% IgM/IgG positivity suggestive of etiologic association. Prospective, comprehensive studies are needed to confirm association of these agents, singly or in combination with AES in Gorakhpur region.

Envelope specific T cell responses & cytokine profiles in chikungunya patients hospitalized with different clinical presentations.

BACKGROUND & OBJECTIVES: Since the 2006 massive outbreaks, chikungunya (CHIK) is a major public health concern in India. The aim of this study was to assess envelope specific immune responses in patients with chikungunya infection. METHODS: This study included 46 hospitalized patients with chikungunya virus infection (encephalitis, n=22, other systemic involvement, OSI, n=12, classical, n=12) and six controls from Ahmedabad city, Gujarat, India. T cell responses and the levels of Th1, pro/ anti-inflammatory cytokines against the CHIK virus envelope antigens were assessed by lymphocyte proliferation assay and by cytometric bead array in flow cytometry, respectively. RESULTS: Lymphoproliferative response was uniform among the patients. Comparisons of cytokines revealed significantly higher levels of interleukin (IL)-4 and IL-5 in encephalitis, OSI and classical patients versus controls. The levels of tumour necrosis factor (TNF)-α were higher in classical patients categories compared to the controls. Interferon (IFN)-γ levels were lower in encephalitis patients versus control. INTERPRETATION & CONCLUSIONS: Our findings showed recognition of T cell epitopes on the envelope region of chikungunya virus by all patient categories. Lower level of IFN-γ may be associated with the severity of disease in these patients.

Effect of pregnancy on anti-HEV antibody titres, plasma cytokines and the corresponding gene expression levels in the PBMCs of patients presenting with self-recovering clinical and subclinical hepatitis E.

High mortality in pregnant women (PR) is a characteristic of hepatitis E in developing countries. To understand the pathogenesis of HEV infection in self-limiting disease during pregnancy, we compared clinical (PR-patients) and subclinical-HEV-infections in pregnant women in the first (SC-PR-1) and later (2nd and 3rd, SC-PR-2+3) trimesters with the respective healthy controls and acute non-PR patients. The SC-PR-2+3 exhibited lower ALT, bilirubin levels, anti-HEV-IgM/IgG titres than the acute-PR/non-PR-patients (p<0.05-0.0001). IFNγ/IL4ratios indicated Th2/Th1 bias in non-PR and PR-patients respectively. Raised levels of 10/20 plasma cytokines in the non-PR-patients reflect predominant inflammatory response, unaltered- IFNγ/reduced-IFNα responses and a robust chemokine secretion. On contrary, the acute-PR-patients exhibited drastic reduction in majority of the cytokines relative to in the non-PR-patients. Importantly, diminished or unaltered response was noted in the acute-PR-group when compared to the corresponding controls. The only exception was sIL2RA, increasing in both patient categories. Of the 14 genes evaluated, the expression of IFNγ/IL10/IL1A/IL7/CCL2/CCL3/CXCL8/CXCL10 was higher in the non-PR patients. Of these, the expression of IFNγ/IL10/IL1A/CCL2/CCL3/CXCL8 and, additionally, IL2/IL6/TNF genes was higher in the clinical-PRs. Almost identical pattern was noted in the control-PR-2+3 category indicating no influence of HEV infection. Comparison of patient-categories identified significant elevation of IFNγ(P<0.001), CCL2(p<0.01), CXCL8(P<0.05), IL1B(p<0.05) and IL10(P<0.0001) and decrease in CXCL10(<0.05) in the PR-patients. The results suggest antibody-dependent disease severity and impaired immune response in the PR patients. Higher expression of cytokine-genes in the PBMCs did not correlate with the plasma-cytokine levels in the PR-patients.

Genomic characterization of Nipah virus, West Bengal, India.

An intrafamilial outbreak in West Bengal, India, involving 5 deaths and person-to-person transmission was attributed to Nipah virus. Full-genome sequence of Nipah virus (18,252 nt) amplified from lung tissue showed 99.2% nt and 99.8% aa identity with the Bangladesh-2004 isolate, suggesting a common source of the virus.

Role of host immune response and viral load in the differential outcome of pandemic H1N1 (2009) influenza virus infection in Indian patients.

BACKGROUND: An unusually high number of severe pneumonia cases with considerable mortality is being observed with the pandemic H1N1 2009 virus infections globally. In India, all mild as well as critically ill cases were admitted and treated in the government hospitals during the initial phase of the pandemic. The present study was undertaken during this early phase of the pandemic. METHODOLOGY: The role of viral load and host factors in the pathogenesis were assessed by examining 26 mild (MP), 15 critically ill patients (CIP) and 20 healthy controls from Pune, India. Sequential blood and lung aspirate samples were collected from CIP. Viral load and cytokines/chemokine levels were determined from the plasma and lung aspirates of the patients. TLR levels were determined by staining and FACS analysis. Gene profiling was done for both cells in the lung aspirates and PBMCs using TaqMan Low Density arrays. Antibody titres and isotyping was done using HA protein based ELISAs. PRINCIPAL FINDINGS: 13/15 critically ill patients expired. All plasma samples were negative for the virus irrespective of the patient's category. Sequential lung samples from CIP showed lower viral loads questioning association of viral replication with the severity. Anti-rpH1N1-09-HA-IgG titres were significantly higher in critically ill patients and both categories circulated exclusively IgG1 isotype. Critically ill patients exhibited increase in TLR-3, 4, 7 and decrease in TLR-2 expressions. The disease severity correlated with increased plasma levels of IL1RA, IL2, IL6, CCL3, CCL4 and IL10. Majority of the immune-function genes were down-regulated in the PBMCs and up-regulated in the cells from lung aspirates of critically ill patients. No distinct pattern differentiating fatal and surviving patients was observed when sequential samples were examined for various parameters. CONCLUSIONS: Disease severity was associated with pronounced impairment of host immune response.

Complete genome sequences of hepatitis C virus subtype 3i and 3a subtype isolates from India.

AIM: Hepatitis C virus (HCV), a major causative agent of chronic hepatitis, is classified into six major genotypes. Genotype 3 HCV infection is more sensitive to interferon therapy. In India, genotype 3, particularly subtype 3a, HCV infections are common. Three novel HCV subtypes i.e., 3g, 3j, and 3i were identified from India based on partial genomic sequences. This report provides full genome sequences of one isolate each of subtypes 3i and 3a. METHODS: Serum samples positive for subtype 3i and 3a HCV RNA based on core region genomic sequences were studied. Complete HCV genomes were amplified as 11 overlapping PCR fragments and sequenced. RESULTS: The complete genomic sequence of Indian HCV 3i isolate clustered with other genotype 3 sequences, and was closer to subtypes 3b and 3a (80.5% and 79.1% [SD 0.4%] nucleotide identity). Nucleotide similarities were the highest in the core region (86.1-88.7%), and the least in the E2 region (69.4-70.7%). Phylogenetic tree analysis confirmed the existence of a separate subtype 3i. The Indian HCV 3a isolate's complete genomic sequences clustered with previously known genotype 3a sequences with a nucleotide similarity of 91.1% (SD 0.2%). Neither isolates showed evidence of recombination of different HCV genotypes. CONCLUSION: The information on complete genomic sequences of the genotype 3 HCV isolates should be helpful in future studies on HCV evolution and classification, and for development of newer therapeutic and preventive strategies against this infection.